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Abstract. 

We analyze the contributions of Nambu-Goldstone modes to the shear viscosity of 
a superfluid atomic Fermi gas close to unitarity. We show that the low temperature 
experimental values of the shear viscosity to entropy ratio, 77/s, are reproduced 
considering an effective shear viscosity which takes into account the collisions among 
Nambu-Goldstone bosons and the finite size of the optical trap. We predict that, for 
T < O.lTp, r]/s should correlate with the size of the trap and linearly decrease with 
decreasing temperature. For T > O.lTp, wc find agreement with the experimental 
data assuming that the Nambu-Goldstone modes have an anomalous dispersion law. 
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The knowledge of the transport coefficients of a fluid may provide, among other 
things, a key to understand the microscopic physics of the system [H |2]. After the 
conjecture [3], derived by string theory methods, of the existence of a universal bound 
of the shear viscosity to entropy ratio rj/s > h/{4:7rkB), special attention has been 
drawn to those fluids that could saturate the bound, as it has been claimed that their 
microscopic degrees of freedom cannot be identified with well-defined quasiparticles. The 
fluid produced in heavy- ion collisions seems to have extremely small shear viscosity [1] , 
close to the universal bound, and almost ideal hydrodynamical flow was observed in a 
dilute Fermi gas of cold atoms at unitarity [1] , i.e. with an infinite two-body scattering 
length. 

Experiments with trapped cold atomic Fermi gases are able to reach the region of 
infinite scattering length tuning the interaction between the fermionic atoms by means 
of a magnetic- field Feshbach resonance [H Ej El [7]. The strength of the interaction 
between atoms depends on the applied magnetic field and can be measured in terms 
of the s-wave scattering length. By varying the magnetic-field controlled interaction, 
fermionic pairing is observed to undergo the Bose-Einstein condensate to Bardeen- 
Cooper-Schrieffer crossover. It is a remarkable aspect of these fermionic systems that 
for any value of the attractive interaction they are superfluid, provided the temperature 
is sufficiently low. 

The experimental setup used to measure the shear viscosity consists of a degenerate 
mixture of an equal number of spin-l/2-up and spin-l/2-down optically trapped 
fermionic atoms. In particular we shall refer to the results reported in Refs. [HI |9] 
with iV ~ 2.0 X 10^ ^Li atoms in a magnetic field of -B = 840 G, close to the s-channel 
Feshbach resonance at -B = 834.15 G [10]. The unperturbed density in the trap is given 
by (in units h = ks = 1, hereafter used) 



where m is the mass of the ^Li atom and Ep is the Fermi energy of N free fermions in 
a harmonic oscillator potential. In the experimental setup of Refs. [HI [9], the transverse 
frequencies of the trap are Ux — 27r x 1778 s~^, ~ 27r x 1617 s~^ and the net axial 
frequency is cu^ ~ 27r x 76 s~^, corresponding to a Fermi temperature Tp ~ 2.4 10~^K. 

Values of the viscosity can be extracted at both high and low temperatures by 
studying either the anisotropic expansion of the atomic cloud or the breathing mode 
damping [TT]. At high T a dimensional analysis, or a much careful study based on 
a Boltzmann equation, predicts the scaling r] oc T'^/^ [12], in agreement with the 
experimental results [11] . At low T, below the superfluid phase transition, the theoretical 
and experimental study of the transport properties of the system is quite challenging. 

Superfluidity is a phenomenon that occurs after the appearance of a quantum 
condensate that spontaneously breaks a global symmetry of the system [13]. First 
discovered in ^He, it was soon realized that the phenomenon occurs in both bosonic 
and fermionic systems. The property of superfluidity follows from the existence of 
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Nambu-Goldstone modes, hereafter referred to as phonons, associated to the breaking 
of the particle number conservation law. In the present manuscript, we employ the 
same mechanism that explains the behavior of rj in ^He (in the regime where phonons 
dominate) to the evaluation of the viscosity in the superfluid phase of the cold Fermi 
atoms in the unitarity limit. 

In our analysis we neglect the contribution of fermionic modes to the shear viscosity. 
This approximation should be reliable for very low temperatures, T < O.lTp, where the 
density of superfluid fermionic modes is exponentially suppressed and the normal fluid 
corona around the center of the trap is small. Actually, the results of [H] show that the 
fermionic contribution to rj is above the experimental values for any value of T < T^. 
Since 7] is directly proportional to the typical time for transport of momentum in the 
direction orthogonal to the flow, r^, there should be modes with a shorter t±. We 
assume that at very low temperature this mode is the phonon. 

Let us briefly review some well known facts about the computation of rj in superfluid 
^He. Landau postulated a Hamiltonian that dictates the dynamics of the phonons [13] 
and at very low temperature Landau and Khalatnikov computed the shear viscosity 
due to binary collisions of phonons, finding the scaling rj oc [15]. Experimental 

measures of rj showed however a different scaling, close to 1/T [16]. Then, Maris had the 
great insight to formulate that the shear viscosity in ^He at low temperature was instead 
dominated by 1 ^ 2 processes (hereafter 3ph) small-angle collisions [17], only possible 
due to an anomalous dispersion law of these collective modes. He further realized that 
the knowledge of the T scaling behavior of r] could provide information on the momentum 
behavior of the phonon dispersion law. A detailed analysis of the experimental results 
has lead to the conclusion that a complex hierarchy of time-scales exists and that in order 
to properly reproduce the experimental results the transition from the hydrodynamic 
to the ballistic regime of phonons should be taken into account [TH [191 1201 l2T] . 

In the superfluid phase of the cold Fermi system at unitarity effective field theory 
techniques have been applied to construct the Lagrangian of the phonons in an energy 
and momentum expansion [22]. At the lowest order in the expansion the Lagrangian 
reads 

£lo = P(^), X = /i-\/(r)-9o</.-^|^, (2) 

where P is the pressure of a superfluid at T = 0, see e.g. [22], is the phase of the 
condensate and V^(r) is the trapping potential. For vanishing trapping potential, the 
formal expression of £lo can be related by a Legendre transformation to the Hamiltonian 
proposed by Landau to derive the phonon self-interactions. The Lagrangian at the next 
to leading order (NLO) in a momentum expansion has as well been constructed in [22] 
and leads to a phonon dispersion law of the form Ep = Csp{l + '~fp'^) where is the speed 
of sound, 

^ = -n'^L+^c^^, (3) 
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where kp is the Fermi momentum, and ^ ^ 0.37 [221 1211 [2S] is the Bertsch number. The 
two parameters Ci and C2, related to the momentum dependence of the static density 
and to the transverse response functions, can be evaluated by the e-expansion technique 
[26] . finding ci ~ —0.021, and C2/C1 = C(e^). They can also be obtained by a fit of 
Monte Carlo numerical simulations [27] , which gives similar results. The interesting fact 
of these results is that 7 > 0, and therefore 3ph small- angle collisions are kinematically 
allowed. The contributions of these processes to the bulk viscosity coefficients and to the 
thermal conductivity of ultracold atomic systems have been computed in Ref. [28] and 
in Ref. [29], respectively. Nonetheless, the results of Ref. [30], obtained by combining 
the e-expansion technique with a gaussian expansion around the mean field, suggest that 
7 < 0. But, considering gaussian fiuctuations close to unitarity might not be enough to 
determine the dispersion law of the collective modes, because the procedure seems to be 
not self-consistent [31]. Therefore, in the present manuscript we refrain from considering 
7 < 0; we shall consider how our results would change in that case in a future work. 

The shear viscosity due to phonon binary collisions in the cold Fermi fiuid has been 
computed, finding a scaling ?74pii oc 1/T^ [321 [33], the same as in ^He, which however is 
not consistent with experimental values of rj. This fact motivates the inclusion of 3ph 
small-angle collisions which are kinematically allowed if the phonon dispersion law is of 
the form 



where il^{p) is positive in a certain range of momenta. We shall assume that il){p) ^ 1, 
and carry out the computation of the shear viscosity in a perturbative expansion in 

For the evaluation of the contribution of 3ph proceses to the shear viscosity, r/sph, 
we consider the phonon transport equation amended by a collision term that takes into 
account 1^2 splitting and joining process. In order to simplify the analysis, in this 
calculation we neglect the effect of the trapping potential, meaning that we evaluate the 
scattering matrix associated to these processes from the Lagrangian in Eq. (|2|) taking 
V^(r) = 0. We obtain that the square of the scattering amplitude reads 



where p is the mass density of the system, u = and p,p', k' refer to the modulus 
of the momenta of the three phonons participating in the process. The expression of 
this scattering amplitude is formally equivalent, as one could naively expect, to the 
corresponding expression obtained for 3ph interactions in ^He. What changes are the 
explicit values of p and Cg in the two different systems. 

We have solved the Boltzmann transport equation using variational methods, and 
find a solution that fully agrees with that proposed to explain the phonon viscosity due 
to small-angle collisions in ^He [3l]. In particular, we find that the contribution of 3ph 



Ep = Csp{l + i^{p)) , 



(4) 




(5) 
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processes to the shear viscosity coefficient is given by 



where 



M=^^i±^ l*dp'dk'{p'k'{p' + k')r\M\' 

W) + h{k') - b{p' + k')f fp^,,{l + fp){l + 47) , (7) 

where /p*^ is the equihbrium Bose-Einstein distribution function for the phonon of 
momentum p, and we have defined b{p) = pip{p)- The star on the top of the integral 
indicates that the integration is constrained by the energy and momentum conservation 
laws. From the expression above one clearly sees that different forms of the phonon 
dispersion law may produce different T scalings of r/sph. 

Assuming that ip{p) = 7P^, where 7 > is given by Eq. ([3D, one finds ?73ph oc 
which is the same scaling found considering binary collisions. This result is 
in disagreement with the experimental data both at low temperature, because ?73ph 
diverges, and at high temperature, because one finds that ?73ph < 1/ (47r) and below the 
experimental values. This last result motivated us to consider, as in the ^He system, a 
dispersion law with 

with 6 > 0, corresponding to phonons with anomalous dispersion. Upon substituting 
Eq. & in Eq. ^ we find that 

where ipmsix = ji and /3 = — while 10) is, up to some normahzation constants, the 
integral that appears in Eq. ([71) after rescaling the momenta so as to make the integral 
over adimensional variables. We then find that in the region, 0.1 Tp ^ T < 0.17Tp, 

~ 2046.4, independent of T, so that ?73ph oc 1/T. For T < 0.1 Tp only the linear 
and cubic terms in momenta in the phonon dispersion law are relevant, and we find that 
?73ph oc 1/T^, in agreement with the result obtained with vanishing 5. 

As we shall see, the shear viscosity in Eq. ([9]) gives a good agreement with the 
experimental data in the range 0.1 Tp ^ T < T^, but it still diverges when the 
temperature is further decreased. However, the measured shear viscosity of the system 
does not diverge and we shall show that this behavior can be reproduced including 
finite size effects, which start to play a relevant role at very low temperature. Indeed, 
we find that at T ~ O.lTp the phonon mean free path, £ph, becomes larger than the 
size of the system. Then, phonons are ballistic and cannot be described by means of 
hydrodynamics. However, it is still possible to define an effective ballistic viscosity, 
which characterizes the damping of the fluid oscillations, in the transition regime. The 
existence of the effective viscosity is linked to the presence of a boundary and to the fact 



|2 



Shear viscosity in a superfluid cold Fermi gas at unitarity 



6 



that quasiparticles can scatter at the boundary. In experiments with ^He it has been 
shown that the effective shear viscosity that takes into account the balhstic phonons 

?7ball = ^PphCsO , (10) 

leads to excellent agreement with the experimental data [TH [191 1201 [21]. Here, 
Pph = ^45^"^ is the density of the normal fluid component and a is the typical length 
scale of the system, which correlates with the size of the device used to measure the 
shear viscosity, see e.g [18]. Analogous equations can be derived for other fluids, see e.g. 
[35]. In the following we shall assume that the same expression for the effective shear 
viscosity of the superfluid phonons given in Eq. ( fTOl) holds for the phonons of the cold 
Fermi gas. 

In the cold Fermi gas the relevant dimension is the shortest dimension of the trap 
and since is the shortest radius, we assume that a < R^. The actual value of a 
depends on the extension of the superfluid region at the trap center and on the properties 
of the boundary between the superfluid region and the corona of unpaired fermions. We 
postpone the evaluation of a and the detailed treatment of the trap to future work and 
consider that the system consists of an homogenous superfluid up to some distance a 
which is independent of the temperature. Then we determine the numerical value of a 
from the experimental data. Since a is related to the extension of the superfluid region, it 
should decrease with increasing temperature. However, in the region T < O.lTp, where 
phonons are ballistic, the extension of the superfluid region is not strongly dependent 
on the temperature, see e.g. [36] . 

Assuming that the 3ph collisions and the interactions with the boundary can be 
treated as independent processes, the effective shear viscosity due to phonons is given 

by m 

^cff = (%ph + ^ball)~^ • (11) 

The effective shear viscosity depends on two unknown quantities, iprna.x and a, which 
we shall use as fitting parameters. In principle, the 4ph processes discussed in [321 [33] 
should be included in the definition of the effective shear viscosity coefficient. However, 
?74ph has only been evaluated using the LO dispersion law and therefore cannot be 
consistently included in the present analysis. 

Assuming that the entropy is also dominated by phonons, we plot in Fig. [Tjthe values 
of rj/s, together with the partial contributions ?73ph/s (dotted red hue) and ?7baii/s (black 
dashed line). We present results obtained with a = 0.3Rx and ipmax = 0.3, meaning that 
5 ~ 0.87^ (taking a = (0.2 — OA)Rx and tjjmax = 0.28 — 0.4 still gives good agreement 
with the experimental data below T^, we shall present a more detailed numerical analysis 
elsewhere). The green vertical dashed line approximately corresponds to the transition 
temperature between the normal phase and the superfluid phase, ^ (0.21 — 0.25)Tp. 
It has to be understood that the various plots of the phonon contributions are meaningful 
only for values of T < Tc, because phonons do not exist in the normal phase. The blue 
solid hne corresponds to the effective shear viscosity, ?7eff, and is obtained by Eq. f jTTj) . 
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The agreement of this result with the experimental data below is remarkable, given 
that our calculation depends in a non-trivial way on only two parameters. Moreover, 
notice that the blue solid line in the range 0.1 Tp < T < Tc is almost independent of 
the temperature, and this behavior is due to the presence of a minimum of rj-^ph/s at 
T ^ 0.2 Tp, which can be traced back to the ipip) correction, given in Eq. ([8]), to the 
phonon dispersion law in Eq. (jlj). 

The fact that in the range 0.1 Tp ^ T < 0.2 Tp the blue solid line is close to the 
experimental values suggests that phonons have an anomalous dispersion law. Indeed, 
assuming that the fermionic processes and the bosonic processes are independent, the 
process with the smaller value of t± (the typical time of transfer of momentum in 
the direction orthogonal to the flow) will dominate. Then, considering that t] oc t±, 
it follows that the experimental values of the shear viscosity represent a lower bound 
for the shear viscosity contribution of each component. If the 6p^ correction, given 
in Eq. (j8]), is not included in the evaluation of rjgg, we find for T ^ 0.2Tp values 
below the experimental ones and below the universal bound l/(47r). Therefore, our 
simplified model suggests that phonons should have an anomalous dispersion law. The 
caveat is that with increasing temperature the reduction of the superfluid region and 
the modification of the matter distribution inside the trap should be properly included 
in the evaluation of the shear viscosity coefficient. For T ^ 0.1 Tp the largest part of 
the system is superfluid [36] and the matter distribution is not strongly changed with 
respect to T = 0, thus we expect the contribution of phonons to be dominant. But, 
with increasing temperature the extension of the corona increases [36] , and the fermionic 
contribution to the shear viscosity should not be negligible, although it is not clear at 
which temperature the fermionic and bosonic contributions should become comparable. 
Moreover, it is not clear whether bosons and fermions can be treated as independent 
degrees of freedom, because some relaxation process between phonons and fermions 
might become important for T > 0.1 Tp. The proper way to take into account both 
degrees of freedom is to solve the Boltzmann equation which takes into account the 
distribution of fermions and phonons in the trap. We leave this complicated analysis to 
future work. 

A remarkable result of our analysis is that t]/s should strongly depend on the 
trap size. If experiments are conducted increasing (reducing) the size of the trap, but 
keeping Ep constant, we predict that rj/s aX low temperatures should increase (decrease). 
The effect should be sizable at any temperature below Tc, and more prominent at 
T < 0.1 Tp, where it is independent of the detailed form of the phonon dispersion 
law. Indeed, at such low temperatures the correction, ip{p)., to the phonon dispersion 
law is completely negligible and ?7efr is dominated by the ballistic effective viscosity, 
because the mean free path of phonons is larger than the trap size. The experimental 
investigation of the ultralow temperature regime would also clarify whether the decrease 
oi Tj/ s with decreasing temperature is due to ballistic phonons or to effects associated 
to noncondensed pairs [11]. Indeed, in Ref. [H] it was found that rj/ s mildly decreases 
with decreasing temperature, while we predict that rj/s oc T. Moreover, the fermionic 
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Figure 1. Plot ojrj/ s in units ofh/ks as a function ofT/Tp. We report all the various 
contributions separately as well as the effective shear viscosity obtained with Eq. 
The interaction of the phonons with the boundary ( dashed black line ) is obtained from 
Eq. assuming a = Q.ZRx', the contribution of the 3-ph process (dotted red line) is 

obtained from Eq. (0j considering ^/^max = 0.3. The experimental values and error bars 
were taken from Ulf . The blue solid line corresponds to the effective shear viscosity 
obtained by Eq. ill]) . The green vertical dashed line approximately corresponds to 
the transition temperature between the normal phase and the superfluid phase 0/. 
The actual critical temperature may be lower, Tc ~ 0.21Tp as reported in \36f . The 
horizontal black dashed line corresponds to the universal bound 0'- 



contribution should not scale with the size of the experimental system, while we find 
that r]/s oc Rx- Whether the effective description of the viscosity proposed in the present 
paper remains valid at such a low temperature or whether some unknown mechanism 
prevents rj/s to vanish with decreasing temperature deserves further investigation. 

It is important to note that the modified dispersion law of phonons, Eq. (jl]), 
infiuences all the thermodynamic properties of the system at low T. These include 
the pressure, the entropy or the specific heat, that would get extra T corrections which 
also depend on the values 7 and 6, if a phonon dispersion law of the form suggested 
in Eq. ([8]) is used. Low thermal corrections to the speed of sound and its attenuation 
factor might also be sensitive to the higher order corrections to the phonon dispersion 
law. It was actually the precise measurement of the specific heat of the low T regime of 
superfiuid ^He [3H] that could confirm the form of the phonon dispersion law necessary 
to explain the value of shear viscosity of that system. It would thus be interesting to 
have precise measurement of different low T thermodynamical variables to confirm the 
form of the phonon dispersion law that we suggested in this manuscript. 
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